Polariton Analysis of a Four-Level Atom Strongly Coupled to a Cavity Mode 
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We present a complete analytical solution for a single four-level atom strongly coupled to a cavity 
field mode and driven by external coherent laser fields. The four-level atomic system consists of 
a three-level subsystem in an E1T configuration, plus an additional atomic level; this system has 
been predicted to exhibit a photon blockade effect. The solution is presented in terms of polaritons. 
An effective Hamiltonian obtained by this procedure is analyzed from the viewpoint of an effective 
two-level system, and the dynamic Stark splitting of dressed states is discussed. The fluorescence 
spectrum of light exiting the cavity mode is analyzed and relevant transitions identified. 

PACS numbers: 42.50.-p, 32.80.-t, 42.65.-k 



I. INTRODUCTION 

The interaction of a single mode of the electromagnetic 
field with a single atom has long been at the forefront of 
interest within the quantum optics community. In this 
context, the Jaynes-Cummings model Q and its exten- 
sions have been the main focus of attention, for several 
reasons. It is the simplest possible model, involving a 
single two-level atom interacting with a quantized field 
mode, and therefore is in many cases exactly solvable. 
It has also proven to be experimentally realizable, thus 
allowing direct comparison between theory and experi- 
ment. This line of research has deepened immensely our 
understanding of fundamental quantum phenomena, and 
continues to do so. 

Experimentally, the field of cavity quantum electrody- 
namics (CQED) |2) has been shown to be very promising 
for further studies of fundamental quantum systems. Re- 
cent advances in mirror manufacturing techniques make 
it possible to build high-finesse microcavities in which the 
coupling strength of an atomic transition to a cavity field 
mode can be an order of magnitude larger than the deco- 
herence rates of the system [9. Furthermore, four spec- 
tacular experiments have recently demonstrated that it is 
possible to trap a single atom within a microscopic cav- 
ity using either an independent atomic trap [§, [5) , or the 
field mode itself, containing not more than one photon 
at a time || f?J . 

Within the framework of CQED, the regime of strong 
atom- field coupling is interesting for many reasons. For 
one, it enables the study of strongly coupled quantum 
systems; in particular, the 'atom-cavity molecule' [§. 
Secondly, it is also a very promising candidate for the 
realization of strong optical nonlinearities Q . For exam- 
ple, an approximation to a x' 3 ' (Kerr) nonlinear opti- 
cal system can be achieved using either a single, strongly 
coupled two- level atom, or an ensemble of weakly coupled 
two-level atoms. However, the large atom-field detuning, 
which minimizes atomic spontaneous emission noise, also 
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minimizes the strength of the nonlinearity. Using the 
Kerr-type nonlinearities produced by a single two-level 
atom in a cavity, conditional quantum dynamics have 
been demonstrated by Brune et al. j|] in the microwave 
regime, and by Turchette et al. |l(J in the optical regime. 

The effect of clcctromagnctically induced transparency 
(EIT) jllj] has been used by Schmidt and Imamoglu 
to devise a scheme involving four-level atoms which pro- 
duces a large Kerr nonlinearity with virtually no noise. 
It has been shown by Imamoglu et al. |l3| that if such 
a strong optical Kerr nonlinearity is implemented in a 
CQED setting, then it is possible to realize photon block- 
ade, in which the atom-cavity system mimics an ideal 
two-level system, and effectively acts as a photon turn- 
stile device for single photons. 

The proposal of Schmidt and Imamoglu jl2| is very 
appealing in its use of EIT to substantially reduce de- 
coherence. To utilize the advantages that a CQED en- 
vironment offers, Rebic et al. ||l4|| proposed a model in 
which a single four-level atom is trapped in a high-finesse 
microcavity. They showed that this system (which we 
call the EIT-Kerr system) can effect a near-ideal Kerr 
optical nonlinearity. In such a strongly coupled system, 
the composite excitations can be labeled as "polaritons" , 
which are defined as mixtures of atom and cavity mode 
excitations. For weak to moderate driving, the EIT- 
Kerr system is well-approximated by a two-state system, 
corresponding to the two lowest lying polariton eigen- 
states. How this behaviour changes with the introduc- 
tion of additional atoms was investigated by Werner and 
Imamoglu [jl5| (see also the work of Greentree et al. jlG| ) . 

Analyzing the single-atom EIT-Kerr system theoreti- 
cally is not in general a straightforward task. In the bad 
cavity regime or the good cavity regime, approximate 
solutions are possible, based on the relative sizes of the 
atom-field coupling constant and the decay rates. In par- 
ticular, it is possible to adiabatically eliminate either the 
cavity or the atomic degrees of freedom, respectively. In 
the strong coupling case, neither of these simplifications 
is possible. The 'atom-field molecule' must be truly re- 
garded as a fundamental entity, which exhibits features 
that cannot be explained in terms of individual proper- 
ties of its constituents. The natural basis for analysis of 
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such a system is the polariton basis. In this paper we per- 
form a polariton analysis of the strongly coupled atom- 
cavity system. Although we concentrate on a particular 
atomic configuration, the underlying method is general 
and could be applied to any strongly coupled system. 

Polariton analysis has been used extensively of late to 
study the dynamics of EIT systems |l7], [l8], , but these 
analyses have concentrated on the semiclassical case of an 
atomic gas driven by laser light; in particular, on the dy- 
namics of 'slow' light. Juzeliunas and Carmichacl j2(J 
have refined the analysis of the corresponding 'slow po- 
laritons', and showed that it is possible to reverse a 
stopped polariton by reversing the control beam. How- 
ever, none of the treatments so far have dealt with the 
coupled atom-cavity system. 

In Section II, we outline the bare model, i.e. the Hamil- 
tonian written in terms of atomic and field operators, 
and explain how dampi ng b y reservoirs enters into the 
formulation. In Section III we diagonalise the interac- 
tion Hamiltonian exactly to find a set of basis states for 
subsequent analysis. In Section IV, the driving term and 
damping terms are expressed in terms of the new basis 
set, and the effective Hamiltonian in the polariton repre- 
sentation is found. In Section [v| we apply our results to 
obtain expressions for the dynamic Stark splitting and 
the spectrum of weak excitations in the effective two- 
level system. In Section [vj we illustrate how to use the 
effective Hamiltonian to identify peaks and linewidths in 
the fluorescence spectrum for the light exiting the cavity 
mode. Finally, conclusions and outlook are presented in 



Section VII 



II. BARE MODEL 

The atomic energy levels are shown in Fig. The atom 
is assumed to be coupled to a single cavity field mode 
and this cavity is driven through one of its mirrors by a 
coherent laser field. The interaction picture Hamiltonian 
describing the system in the rotating wave and electric 
dipole approximations is TL = 7io + Ti-d, where 

Ti Q = hSa 2 2 + ftAo-44 + ihgi (aVi2 - 021a) 

+iHVL c (cr 2 3 - 032) + ihg 2 (aV 3 4 - (743a), ( la ) 



TL d = ih£ P 



(lb) 



Here, tr^ represent atomic raising and lowering operators 
(for i ^ j), and energy level population operators (for 
i = j); (a) is the cavity field creation (annihilation) 
operator. Detunings 6 and A are defined from the rel- 
evant atomic energy levels; g\ y2 are atom-field coupling 
constants for the respective transitions, and Q c is the 
coupling field Rabi frequency. The cavity driving field is 
introduced through the parameter £ p , given by 
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FIG. 1: Atomic energy level scheme. The cavity mode cou- 
ples to transitions |1) — » |2) and |3) — > |4), with respective 
coupling strengths gi and 32 • The transition |2) <-> |3) is cou- 
pled by a classical field of frequency ui c and Rabi frequency 
fi c . Spontaneous emission rates are denoted by 7j . Detunings 
5 and A are defined as positive in the configuration shown. 



In this expression, T is the cavity mirror transmission 
coefficient, k is the cavity decay rate, and V is the power 
output of the driving laser. Damping due to cavity decay 
and spontaneous emission is discussed below. 

Assume that the cavity mode subspace has been trun- 
cated at some finite size A. Together with the four 
atomic levels, these span a Hilbert space of dimension 
4 x A. In the absence of driving (or in the limit where 
term ([It]) becomes negligible), Hamiltonian (jla|) takes a 
block-diagonal form, with A blocks on the main diag- 
onal. Each block represents a manifold of eigenstates 
associated with the appropriate term in the Fock ex- 
pansion. The ground, first and second manifolds have 
been analyzed from the viewpoint of photon blockade 
in Refs. JlJ, [l5|, [r|, where this truncation approach was 
found to be very useful. Addition of the driving term ( |Trj| ) 
significantly complicates the analysis. This term couples 
the different manifolds, and the Hamiltonian matrix loses 
its block-diagonal form. Therefore, it is not practical to 
perform a simple analytical diagonalization of the Hamil- 
tonian (|l|), given the large size of the 4 A by 4 A matrix. 

Dissipation can be added to the model by adding an 
anti-Hermitian term to th«j Hamiltonian (|l|). This term 
results from coupling to reservoir modes, and is obtained 
by tracing the system over these modes. In this approach 
we identify collapse operators, each of them correspond- 
ing to one decay channel |2l[ . In the EIT-Kerr case, there 
are the following four collapse operators 



(3) 



Cl = y/jlCTi2, C*2 = v / 72 cr 32, 
C3 = ^/l3CF^.4:, Ci = yJ~K,a, 



(2) 



where jk denote spontaneous emission rates into each of 
the decay channels, and k denotes the cavity intensity de- 
cay rate. The effective non-Hermitian Hamiltonian takes 
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the form 

4 

H eff = H-iHj2ctC k , (4) 
fc=i 

H being given by (|l]). 

III. DRESSED STATES ANALYSIS 

In this Section we solve the eigenvalue problem exactly 
for the Hamiltonian Tio given by Eq. (Jla|), and obtain a 
basis for further calculations. In a strongly coupled sys- 
tem such as the one under analysis, dressed states Q 
represent the natural basis for analysis, since the system 
under consideration should be viewed as an 'atom-cavity 
molecule', rather than the mere sum of its constituent 
parts (atom + cavity mode in this case). We have al- 
ready remarked in Section II on the complexity of the 
problem of finding the exact (with coherent driving in- 
cluded) dressed states. Alsing et al. J2j| succeeded in 
obtaining the exact solution for the case of a two-level 
atom when the driving field is resonant with the cavity 
mode. They recognized that the eigenstates can be ex- 
pressed as a direct product of field and atomic states, 
where the field states are displaced squeezed states, thus 
simplifying the calculation. The driven EIT-Kerr system 
does not have a solution for the field states with similarly 
convenient properties, so the method of Ref. |Q can not 
be consistently applied. Instead we opt for an alternative 
approach which will be outlined in Section |Tv| . 



|e«) = <4 1) |1 J 1)+,4 1) |0 J 3) 



(6a) 

- u± ) |l,l)+/3i 1) |0,2)+ A1 i 1) |0,3), (6b) 
where the coefficients of the bare states are given by 
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yjl + (ffi/fi c ) 2 + (d'V^) 2 



(7b) 



The energies of these eigenstates are given by E^p — 



K^cav + el 1 '), where 



e (1) - 



(8a) 




+ ,si„ 



Note that J2i=± o ~ ^> reflecting the fact that the 
cavity mode is detuned from a one-photon excitation of 
the atom (see Fig. |l|). 



A. Ground and First Manifold States 

We use the notation \ number of photons in cavity mode, 
atomic energy level ) to denote the bare states. The 
ground state is 

l4 0) ) = |o,i>, (5) 

and has energy E^ = 0. Dressed state j belonging to 

the manifold n is denoted as |e~- ). 

There are three first-manifold states, one of them res- 
onant with the cavity mode, the other two non-resonant 



B. Second and Higher Manifold States 

The second and higher manifold states can be written 
in a generic form, 

i4 w >> = 4 n) M> 2) 

+/4 n) |n-l,3>+4 n) |n-2,4), (9) 

with n > 2 being the manifold label. 

There are four states in each manifold, with energies 
= h{nuj cav + e£ ). The coefficients of these states 

are 
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32 Vn - 1 



= n- 
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ffl(n-l) 
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The exact energies of the four states within a given manifold are found to be, in increasing order, 
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= \ hf T-4 £> ■ 
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C 2 
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2A 

IT 




2B + AC + C 3 /4 
y/C 2 /4-2A/3 + D ' 



where the following abbreviations have been used: 

AS - g\n - gl(n - 1) - fl 2 c , C 



4A 

IT 



A + <f, 



25 + AC + C 3 /4 
^C 2 /4-2A/3 + L> ' 



B = A[ ff 2 n + tt 2 ] + 5g 2 {n~l). / G 2 = (.gi^) 2 n(n - 1) , 
Xi = 2A 3 + 9A(BC - G 2 ) + 27(B 2 - C 2 G 2 ) , X 2 = A 2 + 3BC + 12G 2 , 



X = 



X i + \jX\ - 4Xf , y = X 2 /X , D = (2 1/3 r + 2~ 1/3 X) /3 . 

I 



(10a) 
(10b) 
(10c) 

(lOd) 

(11a) 
(lib) 



(12a) 
(12b) 
(12c) 

(12d) 



Note that for the n-th (n > 2) manifold, X^t=i e i = 
A + <5, which is the two-photon detuning of the atom 
from the cavity resonance. These equations are the exact 
eigenstates of Hamiltonian ([la]), and can be rewritten 
using the polariton operators as 

n, _ h Ji) fi)t (i) , ft (l) (i)t (1) 
rto = " e - P- P- + ' le + P+ P+ 



n=2 j = l 



(13) 



where the polariton operators are defined as p^ = 

\et 1} )(e^ 



Index k in the second row of Eq. (JL3h is 



a dummy index, since P^^P^j 



operators p^\. 



A ; ,e$ n) }(e$ n) |. Polariton 

\ e< o^)( e< o±\ can be written in a rela- 
tively simple form. These expressions, as well as a short 
discussion on the statistical properties of polaritons, are 
given in Appendix [a]. 

The Hamiltonian ( [l3|) is written in a frame rotating 
at the cavity frequency uj cav . The energy level structure 
obtained has a form as shown in Fig. |p dressed states 
with the same detuning relative to the cavity resonance 
will be formally degenerate. There are only two degener- 

Sq "*) and the first manifold 



ate eigenstates: ground state \e 





4 2) ) 


!«?') 




4 2 >) 



FIG. 2: Schematic representation of the dressed states in a 
rotating system: ground, first and second manifold. The two 



states | eg '} and |eQ ±J ) with associated energies e uj = and 
6q = are degenerate. 



J 1 ) 



(o) 



state |eQ ). It can be expected that driving this transi- 
tion coherently yields dynamic Stark splitting. This in- 
deed happens, and will be elaborated upon in Section [v|, 
but next we turn to the problem of how to include driving 
and damping into our polariton model. 



IV. DRIVING AND DAMPING TERMS 

The simplest way to include off-diagonal terms (as as- 
sociated with driving and damping) into the effective 
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Hamiltonian, written in the basis of states calculated in with the polariton operators defined by 



Section III, is to express the atomic and field operators 
in terms of operators describing the transitions between 
these states. Here again, we make a distinction between 
first and higher manifolds, a distinction imposed by the 
very different structure of these two groups of manifolds. 

The first manifold dressed states can be expressed 
in terms of bare states as bi = Midi, where bi = 

(|1, 1), |0, 2),|0,3» T and d x = (| 
The transformation matrix is 

gi/n. 



(- 



\ 



JV_ 

N- 
1 

7V_ 



1 

N u 



N 



N+ 



J 1 ) 



/P-c 



N+ 
1 
N + 



J 



where 



Nn = 



N± = 



s[i + (e ± /n c f + ( 9l /n c y 



(14) 



(15a) 
(15b) 



Similarly, we can write for the higher man- 
ifolds b„ = M„d„, n > 2, where b„ = 



l,2),\n - 1,3), |n - 2,4)) 



(|n, 1), \n 

d n =(\ef\\e^),\e^),\e^))\ a nd 



/ (n)* (n)* (n)* (n)*\ 

' a \ a 2 «3 a\ \ 

^n)* ^n)* p{n)* p (n) 



(n)» 

Mi 



M2 M 3 

(n)* (n)* 

a ^3 



M 4 
4 



and 



(16) 



with the coefficients given by Eqs. ([ic]). These expres- 
sions provide all of the information needed for the sub- 
sequent calculations. 



A. External Driving 

Strongly coupled systems are very sensitive to the num- 
ber of photons. In fact, the most interesting regimes 
include one or a few photons. In the system under inves- 
tigation, the effect of photon blockade occurs when the 
dynamics is limited to the exchange of excitation between 
the ground state and the first manifold. It is therefore 
natural to express the field annihilation operator a in 
terms of the transitions it produces between two adja- 
cent manifolds. In general, deexcitation from manifold n 
to manifold n— 1 occurs via the operator aS n \ expressed 
in terms of the bare states as 



P 



and the full annihilation operator would then be given 
by a sum over all manifolds, a — J2™=i a • For the 
transition from the first excited state to the ground state 
we obtain 



(i.o)„(i) 
Po 



^'°V>\ (18a) 



, j = 0, ±. The effective Rabi frequencies 



(0,1) 



can be calculated from the matrix (|14J) as 



(i,o) _ 



f 



n 



(1,0) 



(siAic) 



2 -^ + (e w )2 



(18b) 
(18c) 



The three terms in the expansion (18a) correspond to 
the three transitions between the ground state and the 
three states excited by a single photon. Each transition 
has an associated effective Rabi frequency ft 



(1,0) 



Note 

that the negative sign in (18c) means the driving of the 
off-resonant states is out of phase with the driving of the 
resonant state (see Fig. || (a)). 

There are twelve possible transitions between the first 
and the second manifolds, driven with the effective Rabi 
frequencies fij^ , where 



f] (2 ' X) 



A> (!)* ( 2 ) 
V2« a) 



+^ (1) */3 7 (2) +MrM 



(1)* (2) 



(19) 



with i — 0, ±; j — 1,... ,4 and it follows from the 
Eq. (||) that (3 { 1] ee (see Fig. | (6)). Since the second 
and subsequent manifolds have four states each, there are 
sixteen transitions between the adjacent manifolds, with 
effective Rabi frequencies of driving 



n 



(n,n — 1) 



'no. 



(n— 1)* (n) 



($ n-1) */3 : 



(») 



+M 



(n—1)* (n) 

M 



+ V~r 



~n (»— 1)* (n) 

2 v) v\ ' 



,(20) 



with n > 2 and i, j = 1, ... ,4 (see Fig. y (c]). The 
coefficients in this expression are given in Eqs (Rw. 



Note that the Rabi frequencies f2 



can also be ob- 



tained from ihn\™' n L> = (e\ n L '\'H<i\e K ™ ) ), with Hd given 



^~ 1] \n d \ef) 

by Eq. (|T|). However, the expansion of the operator a in 
terms of contributions to different transitions, Eq. (p^), 
offers a clearer physical picture of the processes involved 
in the dynamics. The driving Hamiltonian can therefore 
be written in terms of the polariton operators as 



o« = |0, 1)(1, 1| 
aW = V2\l, 1)(2, 1| + |0, 2)(1,2| 
a W = Vn|n-1, l)(n, 1| + V^- 
+\n- 2, 3)(n- 1, 3|) + V^~2\n - 3, 4>(n-2, 4| (17c) 



(17a) 

|0, 3)(1, 3| (17b) 
Tfln-2, 2}{n-l, 2\ 



Hd = ih£p (a — a t ) 



=±,0 



»=±,0 j=l 

oo 4 



(2)t 



(«) _ Jn)t 



(21) 



n—3 



6 



(a) 



1 st manifold 



2 nd manifold 



(c) 

Manifold n 



V V V 



V V V 



Ground State 



I s * manifold 



Manifold (n - 1) 



FIG. 3: Transition between the polaritons in the adjacent manifolds. The cavity resonance is located at the center of each 
manifold. Figures represent: (a) Transitions between the ground state and first manifold states; (6) Transitions between first 
and second manifolds states; (c) Transitions between polaritons in manifolds (n — 1) and n for n > 3. 



Expression ( ^l| ) is the expansion of the driving Hamiltonian in terms of the transitions that are permitted to occur 
between any two dressed states, as shown in Fig. |^. The obvious advantage of this expansion over the original form 
of driving is in the strong coupling/low photon number regime. In this regime, the expansion ( pll) can be truncated 
at the order justified by the problem, while still retaining all (but not more!) of the relevant contributions from the 
external coherent driving. We will illustrate this assertion i n Section [v} 



B. Damping by Reservoir Modes 



The remaining part of the dynamics to be expressed 
in the polariton representation is damping by the reser- 
voir modes. In Section II it was explained how the 
damping enters into the effective Hamiltonian; in partic- 
ular after a trace has been performed over the reservoir 
variables. The resulting Hamiltonian operator is anti- 
Hermitian and has the form 



the contributions from the individual manifolds: 



(n)t„(n) 



a) a = ya y "' ' a 

71=1 

= |1,1>(1,1| 



(23a) 



+ £>|n, l){n, l| + (n-l) (\n-l, 2)(n-l, 2| 

n=2 

+ |n-l, 3)(n-l, 3|) + (n-2)|n-2, 4)(n-2, 4|] , 



Air 



-ihna^a — ih("/i + 72)0-22 — ih^a^ . (22) 



<t 22 = X>&> = X>-1, 2)(n-l, 2|, (23b) 

71 — 1 71=1 
OO OO 

°m = E a $ = E \ n - 2 ' 4 ) < n - 2 > 4 i ■ ( 23c ) 



We follow the reasoning of the previous Section and ex- 
pand the relevant field and atomic operators in terms of 



The operator TL res clearly takes a block-diagonal form in 
the dressed state representation, as the operator expan- 



sion (23) includes terms containing every possible dressed 
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level within a given manifold. The diagonal terms corre- 
spond to damping of the dressed states due to their de- 
cay straight into the reservoir. Off-diagonal terms couple 
two different dressed levels in a given manifold. This cou- 
pling arises due to couplings of both levels to the same 
reservoir level. If each of these diagonal blocks is again 
diagonalized, in the presence of damping we get shifts 
appearing on each level. The complex eigenvalues add 
their real part to the energy shift and the imaginary part 
becomes the damping rate. Energies and damping rates 
calculated in this manner will coincide with the exper- 
imentally observed ones (in the absence of driving). It 
was pointed out by Harris [|24j and Imamoglu |25| that 
these cross terms can be essential in creating destructive 
interference between the transition amplitudes of the ap- 
propriate transitions (see also Li and Xiao [p6|). 

The contribution of the off-diagonal terms to the 
eigenenergies and damping rates (diagonal elements of 
Tires) is very small, and we will ignore their contribution 
to the eigenvalues in the rest of this paper for simplic- 
ity, though we leave them in a general expression for the 
damping Hamiltonian. 

We write the damping Hamiltonian in a form that em- 
phasizes the diagonal and off-diagonal contributions, 



Tt r 



•t rC 1 ) (!) + f 1 ) 
i=±,0 

-*» E ^pfp^ 

i^fe=±,o 

oo 4 

n=2j=l 
oo 

n=2 j^k 



(24) 



where 



T-i(n) (n)* (n) 



(n)* o(n) 



+ [(n-l) K + ll+l2 }f3 n >*p k 



+(n - l)^f>^ ] + [(" - 2)« + 73 ] vfKW (25a) 



in a simple closed form as 



r (L _ 

J- n — 



->(!) 



i + (ffi/n c )^ 

(7i + 72)(d 1} ) ; 



«5i 



5? 



(4")' 



(27a) 
(27b) 



As before, all of the T's could have also been calcu- 



lated from —ihT^ = 



(e^ l '\Ti rea \e£' > ), but again, the outlined procedure offers 
a deeper physical insight. 



» 
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and 



-mr<$ = 



C. Quantum Jumps 

The effective Hamiltonian (Q) describes the time evo- 
lution of the quantum system between successive jumps. 
The effect of quantum jumps is not included, and the 
proper way to include these is the subject of the quan- 
tum trajectories approach Q. Here, we briefly describe 
the transformation of collapses into the dressed state ba- 
sis. 

Quantum jumps are included in the master equation 
for the time evolution of the density matrix p via terms 
of the form CjpCj, where Cj denotes a collapse opera- 
tor from the set (|1) . Each of the collapse operators can 
then be expressed in terms of polariton operators, and 



r5 n) pj;> can 



a new set of collapse operators = 
be obtained. Note that the effective master equation re- 
sulting from the polariton expansion will contain cross 
terms in the collapse operators, giving damping terms of 

the form T\ 



'i] r-'ik ^ik ^ij r l^^ik 

These cross terms have a very important role in modi- 
fying the emission rate from the polariton states to the 
reservoir. Such terms have been studied and well un- 
derstood for the case of the modification of spontaneous 
emission in multilevel atoms 129L 30 . 



V. DYNAMIC STARK SPLITTING 



For each manifold (n), matrix 1^^ is a positive definite 
matrix, so we can write [E7l p8| 



rW 
L jk 

cos 9jk 



cos 9 jk 

Hi ■ fJ'k 

\Pj\\l*k\ 



r (n)p(n) 
1 A j 1 kk 



(26a) 
(26b) 



where fj,j t k can be thought of as the effective dipole mo- 
ments of the transitions contributing the off-diagonal 
term. Furthermore, we note that the diagonal matrix 
elements belonging to the first manifold can be written 



It was proven earlier |l4|, [lj] that the EIT-Kerr sys- 
tem can behave as an effective two-level system, with 



the states |e 



(0)\ 



and 



being its ground and excited 



states. The two-level approximation is best for large val- 
ues of effective dipole coupling, i.e. (g 1 /il c ) 2 ^> 1 and 
g-2 ^> K. The two states of the effective model are cou- 
pled by the external field, with the Rabi frequency of 
the coupling, Oq ' 1 '', given by Eq. ( JlSbj ), and the decay 
rate of the excited state, r^, given by Eq. (27a). It is 



therefore expected that this system exhibits a dynamic 
Stark splitting, characteristic of every driven two-state 
system. We now explore this effect in more detail, using 
the polariton model. 
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FIG. 4: Graphical representation of the dynamic Stark split- 
ting of the effective two-level system. 



Recall that the effective Hamiltonian (Q) is non- 
Hermitian. We truncate the expansion over manifolds of 
the polariton Hamiltonian after the first manifold. More- 
over, we concentrate on the ground and excited state 
of the effective two level system and write its polariton 
Hamiltonian in the reduced form 



H, 



^ = im^( P P- p w)-mr£V^ti ) - (28) 



The eigenvalues of such an effective Hamiltonian are com- 
plex: 



r (i) 



n 



(1,0)2 



(r^/2)" 



(29) 



The real parts of these eigenvalues represent the energies 
of the dressed states, while the imaginary parts represent 
their associated decay rates. We identify two operating 
regimes, depending on the size of fig i i.e. the size of 



• Regime 1: fi 1,0) < T ( 1] /2 



= 0, 



o 
2 



2i^±W(lf'°V2) -n& 



>(1,0) 2 



(30a) 



• Regime 2: n { m > /2 



-.a) 



,(1,0) 2 



r«/2 



(30b) 



The eigenstates corresponding to the Stark-split states 
are 



(IO ± |e 1J » 

i,i)+ 9l /n c \o,3) 



V2 

■±=(\0,1)± 

V2 V V 1 + (SiM) 2 



(31) 



The transition between the two regimes happens at 
°^ = ^0 /2, or J i n terms of the original parameters, 



at 



£ p — 



k/2 



i + ( gi /n c y 



(32) 



Note that there is no contribution of the atomic decay 
rates in this simple model. This absence occurs because 
the excited state in the effective model is a dark state 
with respect to atomic spontaneous emission. The dy- 
namic Stark splitting effect is shown schematically in 
Fig.§. 

Another effective two-level system was predicted by 
Tian and Carmichael J3l| , who studied the case of a two- 
level atom in the cavity, driven on the lower Rabi reso- 
nance. They also predicted Stark splitting |32| compara- 
ble to that presented in this paper. 

For further comparison, the analysis of a two-level 
atom with spontaneous emission rate 7, coupled with 
strength g to the vacuum cavity mode, yields results 
identical to those of Eqs. (p0|), with the replacements 
fig 1,0 '' — > g and r o — > 7 + n. Thus, the dynamic Stark 
splitting found in the EIT-Kerr system can be thought 
of as the exact counterpart to the vacuum Rabi splitting 
characteristic for a two-level atom coupled to the cav- 
ity mode. The distinction of the EIT-Kerr system is in 



the fact that the parameters fig 1 ' ^ and T^' can be ad- 
justed by a simple adjustment of the coupling laser Rabi 
frequency fl c . 

How well can these results describe the full EIT-Kerr 
system, including damping terms, as described by Hamil- 
tonian (§)? Fig. I compares the two solutions. There is 
very good agreement between numerical solution and an- 
alytical approximation, which breaks down only for large 
values of £ p , where truncation after the first manifold is 
not justified any more, since the contribution of states 
from higher manifolds cannot be ignored. 

One additional feature can be seen by looking at the 
eigenenergies in Figure |[ Notice that the Stark splitting 
of eigenenergies does not start at £ p = + but at some 
small, finite value of £ p . This behaviour and the related 
behaviour of the decay rates for weak excitation is shown 
in Fig. H In the limit £ p = 0, the ground and excited 
states of the effective two level system are uncoupled, so 
the decay rates separate accordingly to the decay rates of 
ground state (which is zero) and the excited state T^K 
Increasing the driving strength mixes these two states 
so that their decay rates become approximately equal. 
Once £ p exceeds the value given by Eq. the energy 



.(i) 
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FIG. 5: Numerical (dashed lines) eigenenergies for the full sys- 
tem (cavity mode subspace truncated at 40) compared with 
the anaytical solution of Eq. ( |29[ ) (solid lines). The parame- 
ters are Qj/n = 6, fj/n = 0.1, £I c /k = 2. 



levels shift in opposite directions, giving rise to the Stark 
splitting. For the parameters of Fig. g, this happens at 
£ p = 0.16k. 

In the weak driving regime (Regime 1 above), £ p < 

(k/2)/J 1 + (g 1 /n c ) 2 , and the EIT-Kerr system truly 
behaves as an effective two-level system due to the ab- 
sence of normal mode splitting. This is also the ideal 
photon blockade regime, and will be called a weak driving 

regime. The case £ p > (n/2)/ \J 1 + (gi/£l c ) 2 (Regime 
2) then includes the intermediate and strong driving 
regimes, a study of which will be published elsewhere. 
Once again, we emphasize the similarity of the effective 
two-level system coupled to the driving field and the two- 
level atom coupled to vacuum cavity mode. There is 
an equivalence in the behaviour of the Rabi-split states 
in the latter |33|, |4| to the behaviour of the Stark-split 
states in the former. 

The analysis of the dynamics of the Stark splitting in 
the dressed state basis offers a simple example of the 
convenience of the polariton approach. Once the Hamil- 
tonian is expressed in terms of the polaritons and the re- 
duced effective Hamiltonian is identified, the subsequent 
analysis is considerably simplified. 



VI. FLUORESCENCE SPECTRUM 

The prediction of dynamic Stark splitting in Section |v| 
leads us naturally to an examination of the fluores- 
cence spectrum of the light emitted by the 'atom-cavity 
molecule'. We solve the master equation of the problem 
numerically to obtain the spectrum, and interpret the re- 
sult using the insight provided by the polariton analysis. 




0.1 0.2 

£pl K 



FIG. 6: Eigenenergies and their associated decay rates of the 
Stark split states in the weak excitation regime for the same 
set of parameters as in Fig. |^. 



The master equation of the full atom/cavity system 
may be written in the bare form as 



eff 



i=l 



CipCj 



(33) 



where p is the density matrix of the system, Tl e ff is given 
by Eq. (^) and C, denote the four collapse operators of 
Eq. (^). Using the quantum regression theorem [p5f , we 
solve the master equation and calculate the steady-state 
fluorescence spectrum, 



S F (w) = Re 



lira / drtaUt), a(t + r))e h 



(34) 



Results, for different values of driving £ p , are shown in 
Figs, ffland ||. 

FigTm shows how the central peak (at the frequency of 
driving) splits into a Mollow triplet |36| . The particular 
values of parameters are given in dimensionless units, and 
the small value of n is chosen to produce narrow, highly- 
resolved peaks. For this set of parameters, Eq. (f32|) pre- 
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diets a threshold value for the appearance of Mollow side- 
bands at £ p = 0.0395. The central peak is the result of 

the two transitions IV'-'') an d IV 1 ^) |^+ ) 

between the Stark doublet states in the ground and ex- 
cited states (see Fig. ^). Transitions \4>-^) aim 
^ IV'^ 1 '') cause the appearance of the sidebands at 
frequencies u> = io cav + (e+ — e_) and u> — u> cav + (e_ — e + ) , 
respectively, where e± are given by Eq. ( (30b| ) . 

The lincwidths of these peaks can also be calculated. 
It is straightforward to write the master equation in the 
polariton picture as 

P = -\ (n redP - P Ul ed ) + 2lf ) p( 1) / 4 1)t , (35) 

where Ti. re d is given by Eq. (|2§|). Equations for the den- 
sity matrix elements in the basis spanned by Stark states 
\ip±) of Eq. ( fnj ) can be derived using standard methods 
to give 

r (i) 

P++ = -rWp++ + -|-, (36a) 

P+ - = -(^-i2^y + . 

r (1) 

-±- P -+-rP. (36b) 

From these equations, it is easy to read the spectral 
linewidths of the Mollow spectrum. The central peak will 
have linewidth Tq , while the sidebands have linewidth 
3Iq /2. This is consistent with the results for resonance 
fluorescence p5| . 

Given the complex energy level structure of the atom- 
cavity molecule, it can be expected that transitions other 



than those producing the Mollow spectrum will be seen 
in the fluorescence spectrum. This is indeed true, and 
Fig. H shows the additional sidebands. These peaks are 
relatively small (~ 10~ 3 ), so the associated transitions 
are not expected to contribute significantly to the dy- 
namics. They do, however, cause a departure from the 
ideal two-level behaviour. 

The transitions responsible for the sidebands are iden- 
tified in Fig. ^, where the relevant energy level struc- 
ture is shown. Contributions of the transitions up to 
the third manifold states can be seen. We find peaks at 
the following frequencies: ±Aj w ±2.3, ±A 2 w ±5.7, 
±A3 w ±6.05 and ±A4 w ±6.3. The transitions corre- 
sponding to these peaks can be identified from the en- 

1 A (3) ( 2 ) A ( 3 ) ( 2 ) 

ergy eigenvalues as Ai = e<j — eij , — Ai = e 2 e -> i 

A - J 2 ) 2« A - f (2) =(1) A - 

A 2 - e 3 - e + , -A 2 — e 2 - e_ , A 3 - e 3 - e_ , 

A (2) A (!) ~(0) rn. 

—A3 = e 2 ~~ e + 1 ^4 = £+ £_ • Ihe tiny asymme- 
try between the position of the positive frequency peaks 
and negative frequency peaks arises from the fact that 
the polariton states are asymmetrically detuned from the 
cavity resonance. The source of this asymmetry can be 
traced to the nonzero values of atomic detunings S and A 
(see Fig. |l|). The linewidths of these sidebands can also 
be computed using the same method that produced the 
Eq. fl36|), but the results would be hard to verify, given 
the large degree of overlap between the adjacent peaks, 
seen in Fig. ||. We thus have a complete explanation of 
the fluorescence spectrum. 



VII. CONCLUSION AND OUTLOOK 

In this paper, we have presented an exact solution to 
the eigenvalue problem of the Hamiltonian for a four-level 
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2 nd manifold 



l si manifold 



Ground 



FIG. 9: Schematic depiction of the relevant transitions for 
Fig. [| Full lines denote the energy levels involved, dashed 
lines denote positions of the cavity resonances. A/s are ex- 
plained in the main text. 



compared to the treatment in terms of the bare states. 
For example, if the atom has N a levels (degrees of free- 
dom), and quantum field mode can be safely truncated 
at some number N c , the problem in the bare state basis 
has a dimension of at least N a x N c . In the dressed state 
basis, we can identify which N p dressed states (and asso- 
ciated polaritons) participate in the dynamics, effectively 
reducing the dimension of the problem to N p < N a x N c . 
The treatment of dynamic Stark splitting in Section |v| 
provided a simple but extremely successful example of 
such reduction. In other words, from all of the (infinite 
number of) dimensions of Hilbcrt space, the polariton 
approach lets us pinpoint those few dimensions that are 
predominantly involved in the system dynamics. 

Of course, this approach does not guarantee that the 
reduced problem will be analytically solvable. There are 
some general limits on solvability in the dressed state ba- 
sis. For example, if the number of atomic levels N a > 4, 
the diagonalisation of the interaction Hamiltonian is im- 
possible in principle, except perhaps in some special 
cases, since N a determines the order of the polynomial 
of the eigenvalue problem. Also, the size of the reduced 
problem N p can still be impractically large. While noth- 
ing much can be done about the first problem, for the 
second one, we foresee ways to simplify the involved nu- 
merics. In particular, the coupled amplitudes approach 
and the effective master equation look promising. We are 
pursuing this avenue at the moment, and will publish our 
findings elsewhere. 
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atom strongly coupled to a cavity mode. The regime of 
strong coupling CQED presents a difficult problem for 
analytical calculation, as well as for the understanding 
of the physics involved. We have shown that consistent 
application of the polariton approach can offer a signifi- 
cant insight and even enable a relatively simple analyti- 
cal treatment of the physical problem. In particular, the 
problem of an externally driven atom/quantum field sys- 
tem has been reduced to the problem of composite excita- 
tions, transitions between which are effectively driven by 
classical fields of Rabi frequencies r2^™' n_1 \ which were 
calculated exactly 

The polariton approach can be interpreted as a change 
of basis in Hilbert space. It is obvious that such a change 
can in general simplify the analysis of the problem. The 
reason for this is that in the dressed state basis the num- 
ber of degrees of freedom can be significantly reduced, 



APPENDIX A: POLARITON OPERATORS OF 
THE FIRST MANIFOLD 



In this Appendix, polariton operators for the first man- 
ifold states are given explicitly, and their commutation 
relations discussed. 

From expressions (||) - we deduce the form of the 
polariton operators in terms of atomic and field operators 



(i)t _ a f + (01/^0)031 



Pa 



P^ 



Vi + Cgi/^c) 2 

(ffi/^c) a f + i (e±/Q c ) 021 - 031 

\/i + ( gi /n c ) 2 + (e¥/nn c ) 2 



(Ala) 
(Alb) 



It is a well-known fact that the polaritons are neither 
bosons nor fermions. The commutation relations satisfied 
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by the operators p^p and p^ are 



i + ( fll /n c 



,., ,,,, t - _ (gi/^c) 2 + (eg/gj^i "^ b) 



[j4 1) ,i»£ )t ] 



l + ( 3l /0 c ) +(e^/fiO c ) 



with D21 = <J22 — Ciij £*3i = f"33 — en- Therefore, strong 
coupling of bosons and fermions yields an excitation (or 
a quasiparticle) of mixed statistics. 



We can define two limits in which the polaritons 
become dominated by the contribution of their con- 
stituents, according to the ratio <7i/f2 c . When gi/Q c <C 



(A2a) 2, polariton pjp , and its corresponding eigenstate 



become dominated by their photonic contribution, while 
polaritons p± and corresponding eigenstates become 
dominated by their atomic contribution. So, different 
polaritons become either photon-like or atom-like. The 
opposite situation occurs for gi/Cl c 1. The region 
<7i/fi c ~ 1 is a "no-mans land", where photons and atoms 
contribute comparably to the composition of polaritons. 
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